library(dplyr) ; library(ggplot2) ; library(ggthemes) ; library(tidyr)
## Functions 
kable_cust = function(df_aux,ndec=5,cap=''){
  return(df_aux %>% knitr::kable(digits = ndec,caption = cap) %>% kableExtra::kable_styling())
}

Estudo

Novo relatório com análises para Ana Júlia - Los Hermanos

IRGA

UFSCar

CCSM

Rhizotron

dframe_rhiz = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/Rizhotron.txt")

dframe_rhiz_ufscar = dframe_rhiz %>%
  filter(Intituto=='UFSCar') %>% select(-Manejo,-Tratamentos) %>% 
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')),
         Rhizotron_2022=10*Rhizotron_2022) %>% 
  pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
               values_to = 'Rhizotron') %>% 
  mutate(Safra=gsub('Rhizotron_','',Safra))

dframe_rhiz_ccsm = dframe_rhiz %>%
  filter(Intituto=='CCSM') %>% select(-Tratamento,-Tratamentos) %>%  
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV'))) %>% 
  pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
               values_to = 'Rhizotron') %>% 
  mutate(Safra=gsub('Rhizotron_','',Safra))

UFSCar

dframe_summ_ufscar = dframe_rhiz_ufscar %>%
    select(Safra,RegHidrico,Tratamento,Rhizotron) %>%
    group_by(Safra,RegHidrico,Tratamento) %>% 
    summarise(Mean = mean(Rhizotron,na.rm=T),
              SD = sd(Rhizotron,na.rm=T)) %>% 
  mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))

Graficos

Gráfico De Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ufscar %>% ggplot(
      aes(x=Tratamento,
          fill=RegHidrico,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(~Safra)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Option 2
dframe_summ_ufscar %>% ggplot(
      aes(x=RegHidrico,
          fill=Tratamento,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(~Safra)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Gráfico De Barras com erro
Option 1
dframe_summ_ufscar %>% ggplot(
  aes(x=Tratamento,
      fill=RegHidrico,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=Tratamento,y=Mean,
        color=RegHidrico,
        #linetype=RegHidrico,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(~Safra)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Option 2
dframe_summ_ufscar %>% ggplot(
  aes(x=RegHidrico,
      fill=Tratamento,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=RegHidrico,y=Mean,
        color=Tratamento,
        #linetype=RegHidrico,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(~Safra)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ufscar %>% ggplot()+
    facet_wrap(~Safra,scales='free')+
    geom_point(aes(x=Tratamento,y=Rhizotron,
                          color=RegHidrico,
                          shape=RegHidrico),
               position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(x=Tratamento,y=Mean,
                          color=RegHidrico,
                          shape=RegHidrico),
               data=dframe_summ_ufscar,
               position = position_dodge(width = .75))+
    geom_errorbar(aes(x=Tratamento,y=Mean,
                          color=RegHidrico,
                          linetype=RegHidrico,
                      ymin=Mean-SD,ymax=Mean+SD),
                  data=dframe_summ_ufscar,
                  width = 0.2,alpha=0.9,
                  position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()

Option 2
dframe_rhiz_ufscar %>% ggplot()+
    facet_wrap(~Safra,scales='free')+
    geom_point(aes(x=RegHidrico,y=Rhizotron,
                          color=Tratamento,
                          shape=Tratamento),
               position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(x=RegHidrico,y=Mean,
                          color=Tratamento,
                          shape=Tratamento),
               data=dframe_summ_ufscar,
               position = position_dodge(width = .75))+
    geom_errorbar(aes(x=RegHidrico,y=Mean,
                          color=Tratamento,
                          linetype=Tratamento,
                      ymin=Mean-SD,ymax=Mean+SD),
                  data=dframe_summ_ufscar,
                  width = 0.2,alpha=0.9,
                  position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()

Tabela Resumo
dframe_summ_ufscar %>% select(Safra,RegHidrico,Tratamento,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')
Média e Desvio Padrao - Rhizotron
Safra RegHidrico Tratamento Mean ± SD
2021 S BRF_PF+FD 0.01001 ± 0.00631
2021 S H10+FD 0.0105 ± 0.00623
2021 S H10+H152 0.03834 ± 0.01675
2021 TRF BRF_PF+FD 0.02045 ± 0.02714
2021 TRF H10+FD 0.00918 ± 0.01088
2021 TRF H10+H152 0.01859 ± 0.0209
2021 TRV BRF_PF+FD 0.00808 ± 0.00513
2021 TRV H10+FD 0.01402 ± 0.0136
2021 TRV H10+H152 0.02275 ± 0.00243
2022 S BRF_PF+FD 0.03457 ± 0.00308
2022 S H10+FD 0.03134 ± 0.00449
2022 S H10+H152 0.04866 ± 0.00287
2022 TRF BRF_PF+FD 0.02024 ± 0.00499
2022 TRF H10+FD 0.01622 ± 0.00192
2022 TRF H10+H152 0.02094 ± 0.00723
2022 TRV BRF_PF+FD 0.01795 ± 0.00295
2022 TRV H10+FD 0.02425 ± 0.00266
2022 TRV H10+H152 0.02961 ± 0.00519

Modelo

2021
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group =as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]

summ_tk %>% kable_cust()
RegHidrico Tratamento mean Group
S H10+H152 0.03834 a
TRV H10+H152 0.02275 a
TRF BRF_PF+FD 0.02045 a
TRF H10+H152 0.01859 a
TRV H10+FD 0.01402 a
S H10+FD 0.01050 a
S BRF_PF+FD 0.01001 a
TRF H10+FD 0.00918 a
TRV BRF_PF+FD 0.00808 a 
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
##                       Df   Sum Sq   Mean Sq F value Pr(>F)  
## RegHidrico             2 0.000107 0.0000534   0.257 0.7761  
## Tratamento             2 0.001277 0.0006387   3.077 0.0709 .
## RegHidrico:Tratamento  4 0.000847 0.0002116   1.020 0.4236  
## Residuals             18 0.003736 0.0002076                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.01001 0.00832 1.20307 0.24454
RegHidricoTRF 0.01044 0.01176 0.88738 0.38658
RegHidricoTRV -0.00192 0.01176 -0.16355 0.87191
TratamentoH10+FD 0.00049 0.01176 0.04153 0.96733
TratamentoH10+H152 0.02833 0.01176 2.40859 0.02695
RegHidricoTRF:TratamentoH10+FD -0.01176 0.01664 -0.70673 0.48878
RegHidricoTRV:TratamentoH10+FD 0.00544 0.01664 0.32723 0.74727
RegHidricoTRF:TratamentoH10+H152 -0.03019 0.01664 -1.81452 0.08630
RegHidricoTRV:TratamentoH10+H152 -0.01367 0.01664 -0.82149 0.42211
2022
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group = as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]

summ_tk %>% kable_cust()
RegHidrico Tratamento mean Group
S H10+H152 0.04866 a
S BRF_PF+FD 0.03457 b
S H10+FD 0.03134 bc
TRV H10+H152 0.02961 bcd
TRV H10+FD 0.02425 bcde
TRF H10+H152 0.02094 cde
TRF BRF_PF+FD 0.02024 cde
TRV BRF_PF+FD 0.01795 de
TRF H10+FD 0.01622 e 
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
##                       Df    Sum Sq   Mean Sq F value   Pr(>F)    
## RegHidrico             2 0.0017678 0.0008839  49.258 5.01e-08 ***
## Tratamento             2 0.0004839 0.0002419  13.483 0.000264 ***
## RegHidrico:Tratamento  4 0.0002684 0.0000671   3.739 0.021999 *  
## Residuals             18 0.0003230 0.0000179                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.03457 0.00245 14.13354 0.00000
RegHidricoTRF -0.01433 0.00346 -4.14327 0.00061
RegHidricoTRV -0.01661 0.00346 -4.80345 0.00014
TratamentoH10+FD -0.00323 0.00346 -0.93337 0.36298
TratamentoH10+H152 0.01409 0.00346 4.07497 0.00071
RegHidricoTRF:TratamentoH10+FD -0.00079 0.00489 -0.16097 0.87391
RegHidricoTRV:TratamentoH10+FD 0.00953 0.00489 1.94779 0.06721
RegHidricoTRF:TratamentoH10+H152 -0.01339 0.00489 -2.73656 0.01356
RegHidricoTRV:TratamentoH10+H152 -0.00244 0.00489 -0.49902 0.62381

CCSM

dframe_summ_ccsm = dframe_rhiz_ccsm %>%
    select(Safra,Manejo,Copa,PortaEnxerto,Rhizotron) %>%
    group_by(Safra,Manejo,Copa,PortaEnxerto) %>% 
    summarise(Mean = mean(Rhizotron,na.rm=T),
              SD = sd(Rhizotron,na.rm=T)) %>% 
  mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))

Graficos

Gráfico de Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ccsm %>% ggplot(
      aes(x=PortaEnxerto,
          fill=Copa,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(Safra~Manejo,nrow=1)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Option 2
dframe_summ_ccsm %>% ggplot(
      aes(x=Copa,
          fill=PortaEnxerto,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(Safra~Manejo,nrow=1)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Gráfico de Barras com erro
Option 1
dframe_summ_ccsm %>% ggplot(
  aes(x=PortaEnxerto,
          fill=Copa,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=PortaEnxerto,y=Mean,
        color=Copa,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(Safra~Manejo,nrow=1)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Option 2
dframe_summ_ccsm %>% ggplot(
  aes(x=Copa,
          fill=PortaEnxerto,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=Copa,y=Mean,
        color=PortaEnxerto,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(Safra~Manejo,nrow=1)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ccsm %>% ggplot()+
  facet_wrap(~Safra)+
  geom_point(aes(x=PortaEnxerto,
                 y=Rhizotron,
                 color=Copa,
                 shape=Manejo),
             position = position_dodge(width = .75),alpha=0.5)+
  geom_point(aes(x=PortaEnxerto,
                 y=Mean,
                 color=Copa,
                 shape=Manejo),
             data=dframe_summ_ccsm,
             position = position_dodge(width = .75))+
  geom_errorbar(aes(x=PortaEnxerto,
                 y=Mean,
                 color=Copa,
                 shape=Manejo,
                ymin=Mean-SD,ymax=Mean+SD),
                data=dframe_summ_ccsm,
                width = 0.2,alpha=0.9,
                position = position_dodge(0.75))+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()

Option 2
dframe_rhiz_ccsm %>% ggplot()+
  facet_wrap(~Safra)+
  geom_point(aes(x=Copa,
                 y=Rhizotron,
                 color=PortaEnxerto,
                 shape=Manejo),
             position = position_dodge(width = .75),alpha=0.5)+
  geom_point(aes(x=Copa,
                 y=Mean,
                 color=PortaEnxerto,
                 shape=Manejo),
             data=dframe_summ_ccsm,
             position = position_dodge(width = .75))+
  geom_errorbar(aes(x=Copa,
                 y=Mean,
                 color=PortaEnxerto,
                 shape=Manejo,
                ymin=Mean-SD,ymax=Mean+SD),
                data=dframe_summ_ccsm,
                width = 0.2,alpha=0.9,
                position = position_dodge(0.75))+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()

Tabela Resumo
dframe_summ_ccsm %>% select(Safra,Manejo,PortaEnxerto,Copa,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')
Média e Desvio Padrao - Rhizotron
Safra Manejo PortaEnxerto Copa Mean ± SD
2021 Convencional FD BRF_PF 0.00521 ± 0.00401
2021 Convencional H152 BRF_PF 0.00837 ± 0.0044
2021 Convencional FD IAC10 0.0035 ± 0.00158
2021 Convencional H152 IAC10 0.01423 ± 0.00303
2021 Convencional FD IAC5 0.01404 ± 0.01049
2021 Convencional H152 IAC5 0.00987 ± 0.00401
2021 Ecologico FD BRF_PF 0.0088 ± 0.0014
2021 Ecologico H152 BRF_PF 0.01557 ± 0.00379
2021 Ecologico FD IAC10 0.00549 ± 0.00033
2021 Ecologico H152 IAC10 0.03573 ± 0.01341
2021 Ecologico FD IAC5 0.01022 ± 0.01135
2021 Ecologico H152 IAC5 0.01021 ± 0.00289
2022 Convencional FD BRF_PF 0.01011 ± 0.00569
2022 Convencional H152 BRF_PF 0.02145 ± 0.00038
2022 Convencional FD IAC10 0.00386 ± 0.00174
2022 Convencional H152 IAC10 0.00909 ± 0.00083
2022 Convencional FD IAC5 0.00171 ± 6e-05
2022 Convencional H152 IAC5 0.01492 ± 0.00074
2022 Ecologico FD BRF_PF 0.00425 ± 0.00071
2022 Ecologico H152 BRF_PF 0.02376 ± 0.00298
2022 Ecologico FD IAC10 0.00566 ± 0.00022
2022 Ecologico H152 IAC10 0.00891 ± 0.00142
2022 Ecologico FD IAC5 0.00459 ± 0.00013
2022 Ecologico H152 IAC5 0.01106 ± 0.00222

Modelo

2021
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~Manejo+Copa +PortaEnxerto+
            Manejo*Copa+Manejo*PortaEnxerto +
            Copa*PortaEnxerto+
            Manejo*Copa*PortaEnxerto
          )

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~Manejo+Copa +PortaEnxerto+
            Manejo*Copa+#Manejo*PortaEnxerto +
            Copa*PortaEnxerto+#+
            Manejo*Copa*PortaEnxerto
          )

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(Manejo,Copa,PortaEnxerto) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group =as.data.frame.list(mtcomp$`Manejo:Copa:PortaEnxerto`)[,1]

summ_tk %>% kable_cust()
Manejo Copa PortaEnxerto mean Group
Ecologico IAC10 H152 0.03573 a
Ecologico BRF_PF H152 0.01557 b
Convencional IAC10 H152 0.01423 b
Convencional IAC5 FD 0.01404 b
Ecologico IAC5 FD 0.01022 b
Ecologico IAC5 H152 0.01021 b
Convencional IAC5 H152 0.00987 b
Ecologico BRF_PF FD 0.00880 b
Convencional BRF_PF H152 0.00837 b
Ecologico IAC10 FD 0.00549 b
Convencional BRF_PF FD 0.00521 b
Convencional IAC10 FD 0.00350 b 
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
##                          Df    Sum Sq   Mean Sq F value   Pr(>F)    
## Manejo                    1 0.0002372 0.0002372   5.620 0.026137 *  
## Copa                      2 0.0001736 0.0000868   2.056 0.149916    
## PortaEnxerto              1 0.0005457 0.0005457  12.928 0.001454 ** 
## Manejo:Copa               2 0.0002731 0.0001366   3.236 0.057000 .  
## Copa:PortaEnxerto         2 0.0007999 0.0004000   9.476 0.000926 ***
## Manejo:PortaEnxerto       1 0.0001861 0.0001861   4.409 0.046447 *  
## Manejo:Copa:PortaEnxerto  2 0.0001221 0.0000611   1.447 0.255160    
## Residuals                24 0.0010130 0.0000422                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.00521 0.00375 1.38946 0.17745
ManejoEcologico 0.00359 0.00530 0.67625 0.50535
CopaIAC10 -0.00171 0.00530 -0.32305 0.74946
CopaIAC5 0.00883 0.00530 1.66454 0.10901
PortaEnxertoH152 0.00316 0.00530 0.59530 0.55722
ManejoEcologico:CopaIAC10 -0.00160 0.00750 -0.21266 0.83339
ManejoEcologico:CopaIAC5 -0.00741 0.00750 -0.98735 0.33332
ManejoEcologico:PortaEnxertoH152 0.00362 0.00750 0.48201 0.63416
CopaIAC10:PortaEnxertoH152 0.00757 0.00750 1.00905 0.32302
CopaIAC5:PortaEnxertoH152 -0.00733 0.00750 -0.97650 0.33856
ManejoEcologico:CopaIAC10:PortaEnxertoH152 0.01589 0.01061 1.49799 0.14717
ManejoEcologico:CopaIAC5:PortaEnxertoH152 0.00054 0.01061 0.05116 0.95962
2022
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~
            #Manejo+
            Copa +PortaEnxerto+
            #Manejo*Copa+Manejo*PortaEnxerto +
            Copa*PortaEnxerto#+
            #Manejo*Copa*PortaEnxerto
          )

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~
            Manejo+Copa +PortaEnxerto+
            Manejo*Copa+Manejo*PortaEnxerto +
            Copa*PortaEnxerto+
            Manejo*Copa*PortaEnxerto
          )

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(Copa,PortaEnxerto) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group =as.data.frame.list(mtcomp$`Copa:PortaEnxerto`)[,1]

summ_tk %>% kable_cust()
Copa PortaEnxerto mean Group
BRF_PF H152 0.02261 a
IAC5 H152 0.01299 b
IAC10 H152 0.00900 c
BRF_PF FD 0.00718 cd
IAC10 FD 0.00476 de
IAC5 FD 0.00315 e 
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))

par(mfrow=c(2,2))
plot(fit_aov)

par(mfrow=c(1,1))
summ_aov
##                          Df    Sum Sq   Mean Sq F value   Pr(>F)    
## Manejo                    1 0.0000021 0.0000021   0.480 0.495002    
## Copa                      2 0.0004488 0.0002244  50.641 2.44e-09 ***
## PortaEnxerto              1 0.0008701 0.0008701 196.364 4.73e-13 ***
## Manejo:Copa               2 0.0000099 0.0000050   1.122 0.342089    
## Manejo:PortaEnxerto       1 0.0000001 0.0000001   0.017 0.896199    
## Copa:PortaEnxerto         2 0.0001877 0.0000938  21.174 5.02e-06 ***
## Manejo:Copa:PortaEnxerto  2 0.0000871 0.0000435   9.827 0.000763 ***
## Residuals                24 0.0001063 0.0000044                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.00718 0.00107 6.71937 0.00000
CopaIAC10 -0.00242 0.00151 -1.60200 0.11963
CopaIAC5 -0.00403 0.00151 -2.66480 0.01228
PortaEnxertoH152 0.01543 0.00151 10.20594 0.00000
CopaIAC10:PortaEnxertoH152 -0.01118 0.00214 -5.23293 0.00001
CopaIAC5:PortaEnxertoH152 -0.00559 0.00214 -2.61646 0.01378

Potencial Hidrico

UFSCar

CCSM

Desenvolvimento Vegetativo

dframe_dv = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/DesenvVegetativo.txt")

#dframe_qf_ufscar = dframe_qf %>%
#  filter(Instituto=='UFSCar') %>% select(-Manejo) %>% 
#  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
#
#dframe_qf_ccsm = dframe_qf %>%
#  filter(Instituto=='CCSM') %>% select(-Tratamento) %>%  
#  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))

UFSCar

CCSM

Qualidade do fruto

dframe_qf = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/QualidadeFruto.txt")

dframe_qf_ufscar = dframe_qf %>%
  filter(Instituto=='UFSCar') %>% select(-Manejo) %>% 
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))

dframe_qf_ccsm = dframe_qf %>%
  filter(Instituto=='CCSM') %>% select(-Tratamento) %>%  
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
summ_qf <- function(dframe,variable_name = ''){
  if (variable_name==''){
    stop()
  }
  dframe = dframe %>% select(one_of(variable_name),RegHidrico,Tratamento)
  
  
  colnames(dframe)[1]<-'Y'
  out<- list()
  
  dframe_summ <- dframe %>% group_by(RegHidrico,Tratamento) %>% 
    summarise(Mean = mean(Y,na.rm=T),
              SD = sd(Y,na.rm=T)) %>% 
    mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
  
  out$dframe_summ <- dframe_summ
  
  out$barp <- dframe_summ %>% ggplot(
    aes(x=RegHidrico,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_2 <- dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=RegHidrico,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_error <-  dframe_summ %>% ggplot(
    aes(x=RegHidrico,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=RegHidrico,y=Mean,
          color=Tratamento,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_error_2 <-  dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=RegHidrico,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=Tratamento,y=Mean,
          color=RegHidrico,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=RegHidrico,
      y=Y,
      color=Tratamento,
      shape=Tratamento),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=RegHidrico,
      y=Mean,
      color=Tratamento,
      shape=Tratamento),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=RegHidrico,
      y=Mean,
      color=Tratamento,
      shape=Tratamento,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot_2 <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=Tratamento,
      y=Y,
      color=RegHidrico,
      shape=RegHidrico),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=Tratamento,
      y=Mean,
      color=RegHidrico,
      shape=RegHidrico),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=Tratamento,
      y=Mean,
      color=RegHidrico,
      shape=RegHidrico,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))

  return(out)
}

summ_qf_ccsm <- function(dframe,variable_name = ''){
  if (variable_name==''){
    stop()
  }
  dframe = dframe %>% select(one_of(variable_name),Manejo,Copa,PortaEnxerto)
  
  colnames(dframe)[1]<-'Y'
  out<- list()
  
  dframe = dframe %>% mutate(Tratamento = paste(Copa,'+',PortaEnxerto,sep=''))
  
  dframe_summ <- dframe %>% group_by(Manejo,Tratamento) %>% 
    summarise(Mean = mean(Y,na.rm=T),
              SD = sd(Y,na.rm=T)) %>% 
    mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
  
  
  out$dframe_summ <- dframe_summ
  
  out$barp <- dframe_summ %>% ggplot(
    aes(x=Manejo,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_2 <- dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=Manejo,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_error <-  dframe_summ %>% ggplot(
    aes(x=Manejo,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=Manejo,y=Mean,
          color=Tratamento,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
   out$barp_error_2 <-  dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=Manejo,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=Tratamento,y=Mean,
          color=Manejo,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=Manejo,
      y=Y,
      color=Tratamento,
      shape=Tratamento),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=Manejo,
      y=Mean,
      color=Tratamento,
      shape=Tratamento),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=Manejo,
      y=Mean,
      color=Tratamento,
      shape=Tratamento,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot_2 <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=Tratamento,
      y=Y,
      color=Manejo,
      shape=Manejo),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=Tratamento,
      y=Mean,
      color=Manejo,
      shape=Manejo),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=Tratamento,
      y=Mean,
      color=Manejo,
      shape=Manejo,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))

  return(out)
}

### model funcitons = 

dframe_aux = dframe_qf_ccsm
response='Massa'

fit_qf_ccsm <- function(dframe_aux,response){
  dframe_aux = dframe_aux %>% select(Manejo,PortaEnxerto,Copa,one_of(response))
  out_models = list()
  out_models$models <- list()
  
  ## Full Model
  out_models$models$fit_full          = lm(data=dframe_aux,get(response)~Manejo*PortaEnxerto*Copa)
  out_models$models$fit_full_aov      = aov(data=dframe_aux,get(response)~Manejo*PortaEnxerto*Copa)
  

  ## All 2*2 ints
  out_models$models$fit_noInt3        = lm(data=dframe_aux,get(response) ~ Manejo*Copa +Manejo*PortaEnxerto + PortaEnxerto*Copa)
  out_models$models$fit_noInt3_aov    = aov(data=dframe_aux,get(response)~ Manejo*Copa +Manejo*PortaEnxerto + PortaEnxerto*Copa)

  ## Only PortaEnxerto*Copa
  out_models$models$fit_CopPort       = lm(data=dframe_aux,get(response)~ Manejo + PortaEnxerto*Copa)
  out_models$models$fit_CopPort_aov   = aov(data=dframe_aux,get(response)~ Manejo + PortaEnxerto*Copa)

  ## Just main effects
  out_models$models$fit_mainEffect     = lm(data=dframe_aux,get(response)~ Manejo + PortaEnxerto + Copa)
  out_models$models$fit_mainEffect_aov = aov(data=dframe_aux,get(response)~ Manejo + PortaEnxerto + Copa)

  ### Tukey
  out_models$tukey <- list()
  
  
  
  tk <- TukeyHSD(out_models$models$fit_full_aov)
  mtcomp <- multcompView::multcompLetters4(out_models$models$fit_full_aov,tk)
  out_models$tukey$fit_full_aov  <-   
    dframe_aux %>% group_by(Manejo,PortaEnxerto,Copa) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$`Manejo:PortaEnxerto:Copa`$Letters)
  
  
  tk <- TukeyHSD(out_models$models$fit_noInt3_aov)
  mtcomp <- multcompView::multcompLetters4(out_models$models$fit_noInt3_aov,tk)
  out_models$tukey$fit_noInt3_aov <- list(  
    dframe_aux %>% group_by(Manejo,Copa) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$`Manejo:Copa`$Letters),
    
    dframe_aux %>% group_by(Manejo,PortaEnxerto) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$`Manejo:PortaEnxerto`$Letters),
    
    dframe_aux %>% group_by(Copa,PortaEnxerto) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$`Copa:PortaEnxerto`$Letters)
  )
  
  tk <- TukeyHSD(out_models$models$fit_CopPort_aov)
  mtcomp <- multcompView::multcompLetters4(out_models$models$fit_CopPort_aov,tk)
  out_models$tukey$fit_CopPort_aov <- list(  
    dframe_aux %>% group_by(Manejo,Copa) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$`PortaEnxerto:Copa`$Letters),
    
    dframe_aux %>% group_by(Manejo) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$Manejo$Letters),
    
    dframe_aux %>% group_by(Copa) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$Copa$Letters),
    
    dframe_aux %>% group_by(PortaEnxerto) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$PortaEnxerto$Letters)
  )
  
  tk <- TukeyHSD(out_models$models$fit_mainEffect_aov)
  mtcomp <- multcompView::multcompLetters4(out_models$models$fit_mainEffect_aov,tk)
  out_models$tukey$fit_mainEffect_aov <- list(  
    dframe_aux %>% group_by(Manejo) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$Manejo$Letters),
    
    dframe_aux %>% group_by(Copa) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$Copa$Letters),
    
    dframe_aux %>% group_by(PortaEnxerto) %>% 
    summarise(Mean=mean(get(response))) %>% 
    arrange(-Mean) %>% bind_cols(Group = mtcomp$PortaEnxerto$Letters)
  )
  return(out_models)
}

UFScar

Massa

results <- summ_qf(dframe_qf_ufscar,variable_name = "Massa")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 81.66667 7.09460 81.67 ± 7.09
S FD+H10 47.00000 2.64575 47 ± 2.65
S H14+H10 69.33333 3.51188 69.33 ± 3.51
S H152+H10 69.33333 6.11010 69.33 ± 6.11
S H228+H10 47.66667 0.57735 47.67 ± 0.58
TRF FD+BRS_PF 86.66667 8.14453 86.67 ± 8.14
TRF FD+H10 82.06667 4.10041 82.07 ± 4.1
TRF H14+H10 76.00000 2.64575 76 ± 2.65
TRF H152+H10 79.50000 8.50000 79.5 ± 8.5
TRF H228+H10 83.00000 8.18535 83 ± 8.19
TRV FD+BRS_PF 71.66667 1.15470 71.67 ± 1.15
TRV FD+H10 75.33333 1.15470 75.33 ± 1.15
TRV H14+H10 54.66667 8.50490 54.67 ± 8.5
TRV H152+H10 74.66667 3.51188 74.67 ± 3.51
TRV H228+H10 68.33333 7.02377 68.33 ± 7.02

Altura

results <- summ_qf(dframe_qf_ufscar,variable_name = "Altura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 5.80000 0.20000 5.8 ± 0.2
S FD+H10 4.86667 0.11547 4.87 ± 0.12
S H14+H10 5.53333 0.05774 5.53 ± 0.06
S H152+H10 5.40000 0.20000 5.4 ± 0.2
S H228+H10 5.06667 0.05774 5.07 ± 0.06
TRF FD+BRS_PF 5.86667 0.28868 5.87 ± 0.29
TRF FD+H10 5.83333 0.15275 5.83 ± 0.15
TRF H14+H10 5.60000 0.00000 5.6 ± 0
TRF H152+H10 5.70000 0.10000 5.7 ± 0.1
TRF H228+H10 5.73333 0.25166 5.73 ± 0.25
TRV FD+BRS_PF 5.43333 0.11547 5.43 ± 0.12
TRV FD+H10 5.50000 0.10000 5.5 ± 0.1
TRV H14+H10 5.13333 0.11547 5.13 ± 0.12
TRV H152+H10 5.46667 0.32146 5.47 ± 0.32
TRV H228+H10 5.53333 0.28868 5.53 ± 0.29

Largura

results <- summ_qf(dframe_qf_ufscar,variable_name = "Largura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 5.26667 0.11547 5.27 ± 0.12
S FD+H10 4.36667 0.05774 4.37 ± 0.06
S H14+H10 5.00000 0.10000 5 ± 0.1
S H152+H10 5.03333 0.15275 5.03 ± 0.15
S H228+H10 4.40000 0.10000 4.4 ± 0.1
TRF FD+BRS_PF 5.43333 0.15275 5.43 ± 0.15
TRF FD+H10 5.33333 0.11547 5.33 ± 0.12
TRF H14+H10 5.23333 0.05774 5.23 ± 0.06
TRF H152+H10 5.40000 0.10000 5.4 ± 0.1
TRF H228+H10 5.30000 0.17321 5.3 ± 0.17
TRV FD+BRS_PF 5.10000 0.00000 5.1 ± 0
TRV FD+H10 5.16667 0.05774 5.17 ± 0.06
TRV H14+H10 4.63333 0.30551 4.63 ± 0.31
TRV H152+H10 4.93333 0.28868 4.93 ± 0.29
TRV H228+H10 4.96667 0.15275 4.97 ± 0.15

AL

results <- summ_qf(dframe_qf_ufscar,variable_name = "AL")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 1.10333 0.02082 1.1 ± 0.02
S FD+H10 1.11667 0.02517 1.12 ± 0.03
S H14+H10 1.10667 0.01155 1.11 ± 0.01
S H152+H10 1.07333 0.03055 1.07 ± 0.03
S H228+H10 1.15333 0.04041 1.15 ± 0.04
TRF FD+BRS_PF 1.08333 0.02517 1.08 ± 0.03
TRF FD+H10 1.09333 0.02082 1.09 ± 0.02
TRF H14+H10 1.07333 0.01155 1.07 ± 0.01
TRF H152+H10 1.05667 0.00577 1.06 ± 0.01
TRF H228+H10 1.08333 0.02082 1.08 ± 0.02
TRV FD+BRS_PF 1.06667 0.02309 1.07 ± 0.02
TRV FD+H10 1.06667 0.02309 1.07 ± 0.02
TRV H14+H10 1.11000 0.05000 1.11 ± 0.05
TRV H152+H10 1.11000 0.01000 1.11 ± 0.01
TRV H228+H10 1.11333 0.03055 1.11 ± 0.03

Suco

results <- summ_qf(dframe_qf_ufscar,variable_name = "Suco")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 48.76667 4.05010 48.77 ± 4.05
S FD+H10 37.46667 3.35012 37.47 ± 3.35
S H14+H10 45.63333 1.50444 45.63 ± 1.5
S H152+H10 45.13333 2.58908 45.13 ± 2.59
S H228+H10 39.86667 0.46188 39.87 ± 0.46
TRF FD+BRS_PF 51.83333 2.10792 51.83 ± 2.11
TRF FD+H10 51.20000 2.07846 51.2 ± 2.08
TRF H14+H10 46.46667 1.76163 46.47 ± 1.76
TRF H152+H10 45.00000 2.70000 45 ± 2.7
TRF H228+H10 45.63333 3.13900 45.63 ± 3.14
TRV FD+BRS_PF 45.13333 0.05774 45.13 ± 0.06
TRV FD+H10 49.56667 0.75056 49.57 ± 0.75
TRV H14+H10 38.63333 4.06735 38.63 ± 4.07
TRV H152+H10 41.90000 1.60935 41.9 ± 1.61
TRV H228+H10 43.43333 3.76475 43.43 ± 3.76

Acidez

results <- summ_qf(dframe_qf_ufscar,variable_name = "Acidez")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 5.29067 0.38493 5.29 ± 0.38
S FD+H10 5.42190 0.28282 5.42 ± 0.28
S H14+H10 5.01870 0.39804 5.02 ± 0.4
S H152+H10 4.60677 0.51469 4.61 ± 0.51
S H228+H10 5.06987 0.18971 5.07 ± 0.19
TRF FD+BRS_PF 5.46163 0.26236 5.46 ± 0.26
TRF FD+H10 5.14260 0.43562 5.14 ± 0.44
TRF H14+H10 4.97363 0.08585 4.97 ± 0.09
TRF H152+H10 5.16717 0.15245 5.17 ± 0.15
TRF H228+H10 5.51613 0.34844 5.52 ± 0.35
TRV FD+BRS_PF 5.15400 0.32922 5.15 ± 0.33
TRV FD+H10 4.99920 0.20471 5 ± 0.2
TRV H14+H10 5.00657 0.06020 5.01 ± 0.06
TRV H152+H10 4.79327 0.15489 4.79 ± 0.15
TRV H228+H10 5.50537 0.33448 5.51 ± 0.33

BrixCorr

results <- summ_qf(dframe_qf_ufscar,variable_name = "BrixCorr")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 8.90000 0.17321 8.9 ± 0.17
S FD+H10 8.80000 0.10000 8.8 ± 0.1
S H14+H10 8.83333 0.15275 8.83 ± 0.15
S H152+H10 8.63333 0.05774 8.63 ± 0.06
S H228+H10 9.56667 0.20817 9.57 ± 0.21
TRF FD+BRS_PF 9.13333 0.15275 9.13 ± 0.15
TRF FD+H10 9.16667 0.05774 9.17 ± 0.06
TRF H14+H10 8.93333 0.11547 8.93 ± 0.12
TRF H152+H10 8.90000 0.10000 8.9 ± 0.1
TRF H228+H10 8.90000 0.00000 8.9 ± 0
TRV FD+BRS_PF 8.83333 0.20817 8.83 ± 0.21
TRV FD+H10 9.20000 0.10000 9.2 ± 0.1
TRV H14+H10 8.86667 0.20817 8.87 ± 0.21
TRV H152+H10 8.66667 0.05774 8.67 ± 0.06
TRV H228+H10 8.96667 0.25166 8.97 ± 0.25

Ratio

results <- summ_qf(dframe_qf_ufscar,variable_name = "Ratio")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 1.70000 0.10000 1.7 ± 0.1
S FD+H10 1.60000 0.10000 1.6 ± 0.1
S H14+H10 1.73333 0.11547 1.73 ± 0.12
S H152+H10 1.90000 0.20000 1.9 ± 0.2
S H228+H10 1.90000 0.10000 1.9 ± 0.1
TRF FD+BRS_PF 1.66667 0.05774 1.67 ± 0.06
TRF FD+H10 1.76667 0.15275 1.77 ± 0.15
TRF H14+H10 1.80000 0.00000 1.8 ± 0
TRF H152+H10 1.75000 0.05000 1.75 ± 0.05
TRF H228+H10 1.60000 0.10000 1.6 ± 0.1
TRV FD+BRS_PF 1.73333 0.15275 1.73 ± 0.15
TRV FD+H10 1.86667 0.05774 1.87 ± 0.06
TRV H14+H10 1.76667 0.05774 1.77 ± 0.06
TRV H152+H10 1.80000 0.10000 1.8 ± 0.1
TRV H228+H10 1.63333 0.05774 1.63 ± 0.06

IT

results <- summ_qf(dframe_qf_ufscar,variable_name = "IT")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 1.77000 0.18358 1.77 ± 0.18
S FD+H10 1.34333 0.11015 1.34 ± 0.11
S H14+H10 1.64667 0.03512 1.65 ± 0.04
S H152+H10 1.58667 0.09713 1.59 ± 0.1
S H228+H10 1.55333 0.02309 1.55 ± 0.02
TRF FD+BRS_PF 1.93000 0.05292 1.93 ± 0.05
TRF FD+H10 1.91667 0.05508 1.92 ± 0.06
TRF H14+H10 1.69000 0.08185 1.69 ± 0.08
TRF H152+H10 1.63667 0.08505 1.64 ± 0.09
TRF H228+H10 1.66000 0.11533 1.66 ± 0.12
TRV FD+BRS_PF 1.62667 0.04509 1.63 ± 0.05
TRV FD+H10 1.86000 0.02000 1.86 ± 0.02
TRV H14+H10 1.39667 0.18148 1.4 ± 0.18
TRV H152+H10 1.48000 0.06083 1.48 ± 0.06
TRV H228+H10 1.60000 0.17349 1.6 ± 0.17

VitC

results <- summ_qf(dframe_qf_ufscar,variable_name = "VitC")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 28.54000 0.60622 28.54 ± 0.61
S FD+H10 35.68000 1.14582 35.68 ± 1.15
S H14+H10 32.45667 3.28649 32.46 ± 3.29
S H152+H10 37.78333 1.63102 37.78 ± 1.63
S H228+H10 33.36000 1.15763 33.36 ± 1.16
TRF FD+BRS_PF 33.21000 1.82371 33.21 ± 1.82
TRF FD+H10 32.76000 1.84984 32.76 ± 1.85
TRF H14+H10 32.50667 1.82870 32.51 ± 1.83
TRF H152+H10 34.24333 0.89131 34.24 ± 0.89
TRF H228+H10 34.38333 1.21187 34.38 ± 1.21
TRV FD+BRS_PF 33.45333 1.59080 33.45 ± 1.59
TRV FD+H10 35.30333 0.89891 35.3 ± 0.9
TRV H14+H10 38.17667 1.34500 38.18 ± 1.35
TRV H152+H10 37.89000 1.84789 37.89 ± 1.85
TRV H228+H10 35.40333 2.43097 35.4 ± 2.43

CCSM

Massa

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Massa")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Massa')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 94.33333 11.50362 94.33 ± 11.5
Convencional BRF_PF+H152 89.33333 10.01665 89.33 ± 10.02
Convencional IAC10+FD 90.33333 4.72582 90.33 ± 4.73
Convencional IAC10+H152 86.66667 5.13160 86.67 ± 5.13
Convencional IAC5+FD 88.00000 12.12436 88 ± 12.12
Convencional IAC5+H152 81.66667 2.30940 81.67 ± 2.31
Ecologico BRF_PF+FD 77.33333 3.05505 77.33 ± 3.06
Ecologico BRF_PF+H152 83.66667 0.57735 83.67 ± 0.58
Ecologico IAC10+FD 95.00000 11.13553 95 ± 11.14
Ecologico IAC10+H152 90.33333 2.51661 90.33 ± 2.52
Ecologico IAC5+FD 80.33333 3.21455 80.33 ± 3.21
Ecologico IAC5+H152 88.00000 3.46410 88 ± 3.46
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Ecologico FD IAC10 95.00000 a
Convencional FD BRF_PF 94.33333 a
Convencional FD IAC10 90.33333 a
Ecologico H152 IAC10 90.33333 a
Convencional H152 BRF_PF 89.33333 a
Convencional FD IAC5 88.00000 a
Ecologico H152 IAC5 88.00000 a
Convencional H152 IAC10 86.66667 a
Ecologico H152 BRF_PF 83.66667 a
Convencional H152 IAC5 81.66667 a
Ecologico FD IAC5 80.33333 a
Ecologico FD BRF_PF 77.33333 a 
summary(models$models$fit_full_aov)
##                          Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo                    1   61.4   61.36   1.235 0.2774  
## PortaEnxerto              1    8.0    8.03   0.162 0.6912  
## Copa                      2  237.2  118.58   2.388 0.1133  
## Manejo:PortaEnxerto       1  148.0  148.03   2.980 0.0971 .
## Manejo:Copa               2  377.4  188.69   3.799 0.0369 *
## PortaEnxerto:Copa         2   46.7   23.36   0.470 0.6304  
## Manejo:PortaEnxerto:Copa  2   96.1   48.03   0.967 0.3945  
## Residuals                24 1192.0   49.67                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 94.33333 4.06885 23.18426 0.00000
ManejoEcologico -17.00000 5.75423 -2.95435 0.00691
PortaEnxertoH152 -5.00000 5.75423 -0.86893 0.39349
CopaIAC10 -4.00000 5.75423 -0.69514 0.49364
CopaIAC5 -6.33333 5.75423 -1.10064 0.28197
ManejoEcologico:PortaEnxertoH152 11.33333 8.13770 1.39269 0.17648
ManejoEcologico:CopaIAC10 21.66667 8.13770 2.66250 0.01363
ManejoEcologico:CopaIAC5 9.33333 8.13770 1.14692 0.26272
PortaEnxertoH152:CopaIAC10 1.33333 8.13770 0.16385 0.87122
PortaEnxertoH152:CopaIAC5 -1.33333 8.13770 -0.16385 0.87122
ManejoEcologico:PortaEnxertoH152:CopaIAC10 -12.33333 11.50845 -1.07168 0.29452
ManejoEcologico:PortaEnxertoH152:CopaIAC5 2.66667 11.50845 0.23171 0.81872 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC10 92.66667 a
Convencional BRF_PF 91.83333 a
Convencional IAC10 88.50000 a
Convencional IAC5 84.83333 a
Ecologico IAC5 84.16667 a
Ecologico BRF_PF 80.50000 a
[[2]]
Manejo PortaEnxerto Mean Group
Convencional FD 90.88889 a
Ecologico H152 87.33333 a
Convencional H152 85.88889 a
Ecologico FD 84.22222 a
[[3]]
Copa PortaEnxerto Mean Group
IAC10 FD 92.66667 a
IAC10 H152 88.50000 a
BRF_PF H152 86.50000 a
BRF_PF FD 85.83333 a
IAC5 H152 84.83333 a
IAC5 FD 84.16667 a 
summary(models$models$fit_noInt3_aov)
##                     Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo               1   61.4   61.36   1.239 0.2759  
## Copa                 2  237.2  118.58   2.394 0.1111  
## PortaEnxerto         1    8.0    8.03   0.162 0.6906  
## Manejo:Copa          2  377.4  188.69   3.809 0.0354 *
## Manejo:PortaEnxerto  1  148.0  148.03   2.988 0.0957 .
## Copa:PortaEnxerto    2   46.7   23.36   0.472 0.6293  
## Residuals           26 1288.1   49.54                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 93.52778 3.70962 25.21223 0.00000
ManejoEcologico -15.38889 4.69234 -3.27958 0.00296
CopaIAC10 -0.91667 4.97698 -0.18418 0.85530
CopaIAC5 -7.00000 4.97698 -1.40648 0.17142
PortaEnxertoH152 -3.38889 4.69234 -0.72222 0.47661
ManejoEcologico:CopaIAC10 15.50000 5.74692 2.69710 0.01211
ManejoEcologico:CopaIAC5 10.66667 5.74692 1.85607 0.07481
ManejoEcologico:PortaEnxertoH152 8.11111 4.69234 1.72859 0.09574
CopaIAC10:PortaEnxertoH152 -4.83333 5.74692 -0.84103 0.40800
CopaIAC5:PortaEnxertoH152 0.00000 5.74692 0.00000 1.00000 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC10 92.66667 a
Convencional BRF_PF 91.83333 a
Convencional IAC10 88.50000 a
Convencional IAC5 84.83333 a
Ecologico IAC5 84.16667 a
Ecologico BRF_PF 80.50000 a
[[2]]
Manejo Mean Group
Convencional 88.38889 a
Ecologico 85.77778 a
[[3]]
Copa Mean Group
IAC10 90.58333 a
BRF_PF 86.16667 a
IAC5 84.50000 a
[[4]]
PortaEnxerto Mean Group
FD 87.55556 a
H152 86.61111 a 
summary(models$models$fit_CopPort_aov)
##                   Df Sum Sq Mean Sq F value Pr(>F)
## Manejo             1   61.4   61.36   0.981  0.330
## PortaEnxerto       1    8.0    8.03   0.128  0.723
## Copa               2  237.2  118.58   1.896  0.168
## PortaEnxerto:Copa  2   46.7   23.36   0.374  0.692
## Residuals         29 1813.5   62.53
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 87.13889 3.48702 24.98951 0.00000
ManejoEcologico -2.61111 2.63594 -0.99058 0.33008
PortaEnxertoH152 0.66667 4.56558 0.14602 0.88492
CopaIAC10 6.83333 4.56558 1.49671 0.14527
CopaIAC5 -1.66667 4.56558 -0.36505 0.71772
PortaEnxertoH152:CopaIAC10 -4.83333 6.45670 -0.74858 0.46014
PortaEnxertoH152:CopaIAC5 0.00000 6.45670 0.00000 1.00000 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Convencional 88.38889 a
Ecologico 85.77778 a
[[2]]
Copa Mean Group
IAC10 90.58333 a
BRF_PF 86.16667 a
IAC5 84.50000 a
[[3]]
PortaEnxerto Mean Group
FD 87.55556 a
H152 86.61111 a 
summary(models$models$fit_mainEffect_aov)
##              Df Sum Sq Mean Sq F value Pr(>F)
## Manejo        1   61.4   61.36   1.023  0.320
## PortaEnxerto  1    8.0    8.03   0.134  0.717
## Copa          2  237.2  118.58   1.976  0.156
## Residuals    31 1860.2   60.01
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 87.94444 2.88690 30.46326 0.00000
ManejoEcologico -2.61111 2.58212 -1.01123 0.31974
PortaEnxertoH152 -0.94444 2.58212 -0.36576 0.71703
CopaIAC10 4.41667 3.16244 1.39660 0.17245
CopaIAC5 -1.66667 3.16244 -0.52702 0.60193 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

Altura

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Altura")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Altura')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 6.10000 0.40000 6.1 ± 0.4
Convencional BRF_PF+H152 5.96667 0.25166 5.97 ± 0.25
Convencional IAC10+FD 6.06667 0.15275 6.07 ± 0.15
Convencional IAC10+H152 5.56667 0.11547 5.57 ± 0.12
Convencional IAC5+FD 5.90000 0.36056 5.9 ± 0.36
Convencional IAC5+H152 5.80000 0.10000 5.8 ± 0.1
Ecologico BRF_PF+FD 5.70000 0.10000 5.7 ± 0.1
Ecologico BRF_PF+H152 5.83333 0.05774 5.83 ± 0.06
Ecologico IAC10+FD 6.13333 0.15275 6.13 ± 0.15
Ecologico IAC10+H152 6.03333 0.11547 6.03 ± 0.12
Ecologico IAC5+FD 5.90000 0.10000 5.9 ± 0.1
Ecologico IAC5+H152 6.03333 0.05774 6.03 ± 0.06
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Ecologico FD IAC10 6.13333 a
Convencional FD BRF_PF 6.10000 a
Convencional FD IAC10 6.06667 a
Ecologico H152 IAC10 6.03333 a
Ecologico H152 IAC5 6.03333 a
Convencional H152 BRF_PF 5.96667 a
Ecologico FD IAC5 5.90000 a
Convencional FD IAC5 5.90000 a
Ecologico H152 BRF_PF 5.83333 a
Convencional H152 IAC5 5.80000 a
Ecologico FD BRF_PF 5.70000 a
Convencional H152 IAC10 5.56667 a 
summary(models$models$fit_full_aov)
##                          Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo                    1 0.0136 0.01361   0.353 0.55825   
## PortaEnxerto              1 0.0803 0.08028   2.079 0.16225   
## Copa                      2 0.0172 0.00861   0.223 0.80174   
## Manejo:PortaEnxerto       1 0.2025 0.20250   5.245 0.03110 * 
## Manejo:Copa               2 0.4539 0.22694   5.878 0.00836 **
## PortaEnxerto:Copa         2 0.1906 0.09528   2.468 0.10602   
## Manejo:PortaEnxerto:Copa  2 0.0117 0.00583   0.151 0.86059   
## Residuals                24 0.9267 0.03861                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.10000 0.11345 53.76929 0.00000
ManejoEcologico -0.40000 0.16044 -2.49316 0.01996
PortaEnxertoH152 -0.13333 0.16044 -0.83105 0.41414
CopaIAC10 -0.03333 0.16044 -0.20776 0.83717
CopaIAC5 -0.20000 0.16044 -1.24658 0.22458
ManejoEcologico:PortaEnxertoH152 0.26667 0.22690 1.17529 0.25141
ManejoEcologico:CopaIAC10 0.46667 0.22690 2.05675 0.05074
ManejoEcologico:CopaIAC5 0.40000 0.22690 1.76293 0.09064
PortaEnxertoH152:CopaIAC10 -0.36667 0.22690 -1.61602 0.11916
PortaEnxertoH152:CopaIAC5 0.03333 0.22690 0.14691 0.88443
ManejoEcologico:PortaEnxertoH152:CopaIAC10 0.13333 0.32088 0.41553 0.68145
ManejoEcologico:PortaEnxertoH152:CopaIAC5 -0.03333 0.32088 -0.10388 0.91813 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC10 6.08333 a
Convencional BRF_PF 6.03333 a
Ecologico IAC5 5.96667 a
Convencional IAC5 5.85000 a
Convencional IAC10 5.81667 a
Ecologico BRF_PF 5.76667 a
[[2]]
Manejo PortaEnxerto Mean Group
Convencional FD 6.02222 a
Ecologico H152 5.96667 a
Ecologico FD 5.91111 a
Convencional H152 5.77778 a
[[3]]
Copa PortaEnxerto Mean Group
IAC10 FD 6.10000 a
IAC5 H152 5.91667 a
BRF_PF FD 5.90000 a
BRF_PF H152 5.90000 a
IAC5 FD 5.90000 a
IAC10 H152 5.80000 a 
summary(models$models$fit_noInt3_aov)
##                     Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo               1 0.0136 0.01361   0.377 0.54447   
## Copa                 2 0.0172 0.00861   0.239 0.78943   
## PortaEnxerto         1 0.0803 0.08028   2.224 0.14788   
## Manejo:Copa          2 0.4539 0.22694   6.288 0.00592 **
## Manejo:PortaEnxerto  1 0.2025 0.20250   5.611 0.02556 * 
## Copa:PortaEnxerto    2 0.1906 0.09528   2.640 0.09040 . 
## Residuals           26 0.9383 0.03609                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.10833 0.10012 61.00734 0.00000
ManejoEcologico -0.41667 0.12665 -3.28994 0.00288
CopaIAC10 -0.06667 0.13433 -0.49629 0.62386
CopaIAC5 -0.19167 0.13433 -1.42682 0.16553
PortaEnxertoH152 -0.15000 0.12665 -1.18438 0.24698
ManejoEcologico:CopaIAC10 0.53333 0.15511 3.43837 0.00198
ManejoEcologico:CopaIAC5 0.38333 0.15511 2.47133 0.02033
ManejoEcologico:PortaEnxertoH152 0.30000 0.12665 2.36876 0.02556
CopaIAC10:PortaEnxertoH152 -0.30000 0.15511 -1.93408 0.06406
CopaIAC5:PortaEnxertoH152 0.01667 0.15511 0.10745 0.91526 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC10 6.08333 a
Convencional BRF_PF 6.03333 a
Ecologico IAC5 5.96667 a
Convencional IAC5 5.85000 a
Convencional IAC10 5.81667 a
Ecologico BRF_PF 5.76667 a
[[2]]
Manejo Mean Group
Ecologico 5.93889 a
Convencional 5.90000 a
[[3]]
Copa Mean Group
IAC10 5.95000 a
IAC5 5.90833 a
BRF_PF 5.90000 a
[[4]]
PortaEnxerto Mean Group
FD 5.96667 a
H152 5.87222 a 
summary(models$models$fit_CopPort_aov)
##                   Df Sum Sq Mean Sq F value Pr(>F)
## Manejo             1 0.0136 0.01361   0.248  0.623
## PortaEnxerto       1 0.0803 0.08028   1.460  0.237
## Copa               2 0.0172 0.00861   0.157  0.856
## PortaEnxerto:Copa  2 0.1906 0.09528   1.733  0.195
## Residuals         29 1.5947 0.05499
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.88056 0.10340 56.86919 0.00000
ManejoEcologico 0.03889 0.07817 0.49751 0.62258
PortaEnxertoH152 0.00000 0.13539 0.00000 1.00000
CopaIAC10 0.20000 0.13539 1.47723 0.15039
CopaIAC5 0.00000 0.13539 0.00000 1.00000
PortaEnxertoH152:CopaIAC10 -0.30000 0.19147 -1.56684 0.12800
PortaEnxertoH152:CopaIAC5 0.01667 0.19147 0.08705 0.93123 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Ecologico 5.93889 a
Convencional 5.90000 a
[[2]]
Copa Mean Group
IAC10 5.95000 a
IAC5 5.90833 a
BRF_PF 5.90000 a
[[3]]
PortaEnxerto Mean Group
FD 5.96667 a
H152 5.87222 a 
summary(models$models$fit_mainEffect_aov)
##              Df Sum Sq Mean Sq F value Pr(>F)
## Manejo        1 0.0136 0.01361   0.236  0.630
## PortaEnxerto  1 0.0803 0.08028   1.394  0.247
## Copa          2 0.0172 0.00861   0.150  0.862
## Residuals    31 1.7853 0.05759
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.92778 0.08943 66.28055 0.00000
ManejoEcologico 0.03889 0.07999 0.48615 0.63028
PortaEnxertoH152 -0.09444 0.07999 -1.18066 0.24672
CopaIAC10 0.05000 0.09797 0.51036 0.61342
CopaIAC5 0.00833 0.09797 0.08506 0.93276 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

AL

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "AL")
models <- fit_qf_ccsm(dframe_qf_ccsm,'AL')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 1.09000 0.01732 1.09 ± 0.02
Convencional BRF_PF+H152 1.10667 0.02082 1.11 ± 0.02
Convencional IAC10+FD 1.12333 0.01155 1.12 ± 0.01
Convencional IAC10+H152 1.09333 0.01155 1.09 ± 0.01
Convencional IAC5+FD 1.08667 0.02082 1.09 ± 0.02
Convencional IAC5+H152 1.10000 0.01000 1.1 ± 0.01
Ecologico BRF_PF+FD 1.10667 0.03055 1.11 ± 0.03
Ecologico BRF_PF+H152 1.10667 0.01528 1.11 ± 0.02
Ecologico IAC10+FD 1.10333 0.01155 1.1 ± 0.01
Ecologico IAC10+H152 1.15333 0.03786 1.15 ± 0.04
Ecologico IAC5+FD 1.11333 0.02082 1.11 ± 0.02
Ecologico IAC5+H152 1.11000 0.02000 1.11 ± 0.02
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Ecologico H152 IAC10 1.15333 a
Convencional FD IAC10 1.12333 ab
Ecologico FD IAC5 1.11333 ab
Ecologico H152 IAC5 1.11000 ab
Convencional H152 BRF_PF 1.10667 ab
Ecologico FD BRF_PF 1.10667 ab
Ecologico H152 BRF_PF 1.10667 ab
Ecologico FD IAC10 1.10333 ab
Convencional H152 IAC5 1.10000 ab
Convencional H152 IAC10 1.09333 ab
Convencional FD BRF_PF 1.09000 b
Convencional FD IAC5 1.08667 b 
summary(models$models$fit_full_aov)
##                          Df   Sum Sq   Mean Sq F value Pr(>F)  
## Manejo                    1 0.002178 0.0021778   5.124 0.0329 *
## PortaEnxerto              1 0.000544 0.0005444   1.281 0.2689  
## Copa                      2 0.002006 0.0010028   2.359 0.1160  
## Manejo:PortaEnxerto       1 0.000544 0.0005444   1.281 0.2689  
## Manejo:Copa               2 0.000239 0.0001194   0.281 0.7574  
## PortaEnxerto:Copa         2 0.000039 0.0000194   0.046 0.9554  
## Manejo:PortaEnxerto:Copa  2 0.004672 0.0023361   5.497 0.0108 *
## Residuals                24 0.010200 0.0004250                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.09000 0.01190 91.57832 0.00000
ManejoEcologico 0.01667 0.01683 0.99015 0.33198
PortaEnxertoH152 0.01667 0.01683 0.99015 0.33198
CopaIAC10 0.03333 0.01683 1.98030 0.05925
CopaIAC5 -0.00333 0.01683 -0.19803 0.84469
ManejoEcologico:PortaEnxertoH152 -0.01667 0.02380 -0.70014 0.49057
ManejoEcologico:CopaIAC10 -0.03667 0.02380 -1.54031 0.13657
ManejoEcologico:CopaIAC5 0.01000 0.02380 0.42008 0.67816
PortaEnxertoH152:CopaIAC10 -0.04667 0.02380 -1.96039 0.06166
PortaEnxertoH152:CopaIAC5 -0.00333 0.02380 -0.14003 0.88981
ManejoEcologico:PortaEnxertoH152:CopaIAC10 0.09667 0.03367 2.87143 0.00841
ManejoEcologico:PortaEnxertoH152:CopaIAC5 0.00000 0.03367 0.00000 1.00000 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC10 1.12833 a
Ecologico IAC5 1.11167 a
Convencional IAC10 1.10833 a
Ecologico BRF_PF 1.10667 a
Convencional BRF_PF 1.09833 a
Convencional IAC5 1.09333 a
[[2]]
Manejo PortaEnxerto Mean Group
Ecologico H152 1.12333 a
Ecologico FD 1.10778 a
Convencional FD 1.10000 a
Convencional H152 1.10000 a
[[3]]
Copa PortaEnxerto Mean Group
IAC10 H152 1.12333 a
IAC10 FD 1.11333 a
BRF_PF H152 1.10667 a
IAC5 H152 1.10500 a
IAC5 FD 1.10000 a
BRF_PF FD 1.09833 a 
summary(models$models$fit_noInt3_aov)
##                     Df   Sum Sq   Mean Sq F value Pr(>F)  
## Manejo               1 0.002178 0.0021778   3.807 0.0619 .
## Copa                 2 0.002006 0.0010028   1.753 0.1931  
## PortaEnxerto         1 0.000544 0.0005444   0.952 0.3382  
## Manejo:Copa          2 0.000239 0.0001194   0.209 0.8129  
## Manejo:PortaEnxerto  1 0.000544 0.0005444   0.952 0.3382  
## Copa:PortaEnxerto    2 0.000039 0.0000194   0.034 0.9666  
## Residuals           26 0.014872 0.0005720                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.09806 0.01261 87.11126 0.00000
ManejoEcologico 0.00056 0.01594 0.03484 0.97247
CopaIAC10 0.00917 0.01691 0.54203 0.59241
CopaIAC5 -0.00333 0.01691 -0.19710 0.84528
PortaEnxertoH152 0.00056 0.01594 0.03484 0.97247
ManejoEcologico:CopaIAC10 0.01167 0.01953 0.59744 0.55538
ManejoEcologico:CopaIAC5 0.01000 0.01953 0.51209 0.61291
ManejoEcologico:PortaEnxertoH152 0.01556 0.01594 0.97561 0.33825
CopaIAC10:PortaEnxertoH152 0.00167 0.01953 0.08535 0.93264
CopaIAC5:PortaEnxertoH152 -0.00333 0.01953 -0.17070 0.86578 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC10 1.12833 a
Ecologico IAC5 1.11167 a
Convencional IAC10 1.10833 a
Ecologico BRF_PF 1.10667 a
Convencional BRF_PF 1.09833 a
Convencional IAC5 1.09333 a
[[2]]
Manejo Mean Group
Ecologico 1.11556 a
Convencional 1.10000 a
[[3]]
Copa Mean Group
IAC10 1.11833 a
BRF_PF 1.10250 a
IAC5 1.10250 a
[[4]]
PortaEnxerto Mean Group
H152 1.11167 a
FD 1.10389 a 
summary(models$models$fit_CopPort_aov)
##                   Df   Sum Sq   Mean Sq F value Pr(>F)  
## Manejo             1 0.002178 0.0021778   4.034  0.054 .
## PortaEnxerto       1 0.000544 0.0005444   1.009  0.324  
## Copa               2 0.002006 0.0010028   1.858  0.174  
## PortaEnxerto:Copa  2 0.000039 0.0000194   0.036  0.965  
## Residuals         29 0.015656 0.0005398                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.09056 0.01025 106.44243 0.00000
ManejoEcologico 0.01556 0.00774 2.00850 0.05398
PortaEnxertoH152 0.00833 0.01341 0.62122 0.53931
CopaIAC10 0.01500 0.01341 1.11819 0.27266
CopaIAC5 0.00167 0.01341 0.12424 0.90198
PortaEnxertoH152:CopaIAC10 0.00167 0.01897 0.08785 0.93060
PortaEnxertoH152:CopaIAC5 -0.00333 0.01897 -0.17571 0.86175 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Ecologico 1.11556 a
Convencional 1.10000 b
[[2]]
Copa Mean Group
IAC10 1.11833 a
BRF_PF 1.10250 a
IAC5 1.10250 a
[[3]]
PortaEnxerto Mean Group
H152 1.11167 a
FD 1.10389 a 
summary(models$models$fit_mainEffect_aov)
##              Df   Sum Sq   Mean Sq F value Pr(>F)  
## Manejo        1 0.002178 0.0021778   4.302 0.0465 *
## PortaEnxerto  1 0.000544 0.0005444   1.075 0.3078  
## Copa          2 0.002006 0.0010028   1.981 0.1551  
## Residuals    31 0.015694 0.0005063                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.09083 0.00839 130.08659 0.00000
ManejoEcologico 0.01556 0.00750 2.07403 0.04647
PortaEnxertoH152 0.00778 0.00750 1.03701 0.30775
CopaIAC10 0.01583 0.00919 1.72368 0.09473
CopaIAC5 0.00000 0.00919 0.00000 1.00000 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

Acidez

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Acidez")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Acidez')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 5.79277 0.21457 5.79 ± 0.21
Convencional BRF_PF+H152 4.73127 1.12089 4.73 ± 1.12
Convencional IAC10+FD 5.78470 0.14056 5.78 ± 0.14
Convencional IAC10+H152 4.64587 0.14255 4.65 ± 0.14
Convencional IAC5+FD 5.50670 0.55673 5.51 ± 0.56
Convencional IAC5+H152 5.28863 0.23447 5.29 ± 0.23
Ecologico BRF_PF+FD 5.43200 0.25991 5.43 ± 0.26
Ecologico BRF_PF+H152 5.71203 0.41669 5.71 ± 0.42
Ecologico IAC10+FD 5.43063 0.73447 5.43 ± 0.73
Ecologico IAC10+H152 5.57807 0.10602 5.58 ± 0.11
Ecologico IAC5+FD 5.59083 0.15735 5.59 ± 0.16
Ecologico IAC5+H152 5.64403 0.08699 5.64 ± 0.09
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Convencional FD BRF_PF 5.79277 a
Convencional FD IAC10 5.78470 a
Ecologico H152 BRF_PF 5.71203 a
Ecologico H152 IAC5 5.64403 a
Ecologico FD IAC5 5.59083 a
Ecologico H152 IAC10 5.57807 a
Convencional FD IAC5 5.50670 a
Ecologico FD BRF_PF 5.43200 a
Ecologico FD IAC10 5.43063 a
Convencional H152 IAC5 5.28863 a
Convencional H152 BRF_PF 4.73127 a
Convencional H152 IAC10 4.64587 a 
summary(models$models$fit_full_aov)
##                          Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo                    1  0.670  0.6705   3.178 0.08728 . 
## PortaEnxerto              1  0.939  0.9387   4.449 0.04553 * 
## Copa                      2  0.133  0.0666   0.316 0.73232   
## Manejo:PortaEnxerto       1  2.101  2.1011   9.960 0.00427 **
## Manejo:Copa               2  0.013  0.0067   0.032 0.96882   
## PortaEnxerto:Copa         2  0.277  0.1384   0.656 0.52790   
## Manejo:PortaEnxerto:Copa  2  0.545  0.2723   1.291 0.29344   
## Residuals                24  5.063  0.2110                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.79277 0.26518 21.84430 0.00000
ManejoEcologico -0.36077 0.37503 -0.96197 0.34566
PortaEnxertoH152 -1.06150 0.37503 -2.83046 0.00925
CopaIAC10 -0.00807 0.37503 -0.02151 0.98302
CopaIAC5 -0.28607 0.37503 -0.76279 0.45302
ManejoEcologico:PortaEnxertoH152 1.34153 0.53037 2.52944 0.01840
ManejoEcologico:CopaIAC10 0.00670 0.53037 0.01263 0.99003
ManejoEcologico:CopaIAC5 0.44490 0.53037 0.83885 0.40983
PortaEnxertoH152:CopaIAC10 -0.07733 0.53037 -0.14581 0.88529
PortaEnxertoH152:CopaIAC5 0.84343 0.53037 1.59028 0.12486
ManejoEcologico:PortaEnxertoH152:CopaIAC10 -0.05527 0.75005 -0.07368 0.94187
ManejoEcologico:PortaEnxertoH152:CopaIAC5 -1.07027 0.75005 -1.42692 0.16648 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC5 5.61743 a
Ecologico BRF_PF 5.57202 a
Ecologico IAC10 5.50435 a
Convencional IAC5 5.39767 a
Convencional BRF_PF 5.26202 a
Convencional IAC10 5.21528 a
[[2]]
Manejo PortaEnxerto Mean Group
Convencional FD 5.69472 a
Ecologico H152 5.64471 a
Ecologico FD 5.48449 ab
Convencional H152 4.88859 b
[[3]]
Copa PortaEnxerto Mean Group
BRF_PF FD 5.61238 a
IAC10 FD 5.60767 a
IAC5 FD 5.54877 a
IAC5 H152 5.46633 a
BRF_PF H152 5.22165 a
IAC10 H152 5.11197 a 
summary(models$models$fit_noInt3_aov)
##                     Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo               1  0.670  0.6705   3.109 0.08963 . 
## Copa                 2  0.133  0.0666   0.309 0.73705   
## PortaEnxerto         1  0.939  0.9387   4.352 0.04692 * 
## Manejo:Copa          2  0.013  0.0067   0.031 0.96949   
## Manejo:PortaEnxerto  1  2.101  2.1011   9.742 0.00438 **
## Copa:PortaEnxerto    2  0.277  0.1384   0.642 0.53448   
## Residuals           26  5.608  0.2157                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.69897 0.24477 23.28270 0.00000
ManejoEcologico -0.17318 0.30962 -0.55933 0.58072
CopaIAC10 0.00575 0.32840 0.01751 0.98616
CopaIAC5 -0.01850 0.32840 -0.05633 0.95551
PortaEnxertoH152 -0.87391 0.30962 -2.82257 0.00901
ManejoEcologico:CopaIAC10 -0.02093 0.37920 -0.05520 0.95640
ManejoEcologico:CopaIAC5 -0.09023 0.37920 -0.23796 0.81378
ManejoEcologico:PortaEnxertoH152 0.96636 0.30962 3.12114 0.00438
CopaIAC10:PortaEnxertoH152 -0.10497 0.37920 -0.27681 0.78411
CopaIAC5:PortaEnxertoH152 0.30830 0.37920 0.81303 0.42359 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico IAC5 5.61743 a
Ecologico BRF_PF 5.57202 a
Ecologico IAC10 5.50435 a
Convencional IAC5 5.39767 a
Convencional BRF_PF 5.26202 a
Convencional IAC10 5.21528 a
[[2]]
Manejo Mean Group
Ecologico 5.56460 a
Convencional 5.29166 a
[[3]]
Copa Mean Group
IAC5 5.50755 a
BRF_PF 5.41702 a
IAC10 5.35982 a
[[4]]
PortaEnxerto Mean Group
FD 5.58961 a
H152 5.26665 a 
summary(models$models$fit_CopPort_aov)
##                   Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo             1  0.670  0.6705   2.518 0.1234  
## PortaEnxerto       1  0.939  0.9387   3.525 0.0705 .
## Copa               2  0.133  0.0666   0.250 0.7804  
## PortaEnxerto:Copa  2  0.277  0.1384   0.520 0.6001  
## Residuals         29  7.722  0.2663                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.47591 0.22755 24.06467 0.00000
ManejoEcologico 0.27294 0.17201 1.58678 0.12341
PortaEnxertoH152 -0.39073 0.29793 -1.31148 0.19998
CopaIAC10 -0.00472 0.29793 -0.01583 0.98748
CopaIAC5 -0.06362 0.29793 -0.21353 0.83241
PortaEnxertoH152:CopaIAC10 -0.10497 0.42134 -0.24913 0.80502
PortaEnxertoH152:CopaIAC5 0.30830 0.42134 0.73171 0.47022 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Ecologico 5.56460 a
Convencional 5.29166 a
[[2]]
Copa Mean Group
IAC5 5.50755 a
BRF_PF 5.41702 a
IAC10 5.35982 a
[[3]]
PortaEnxerto Mean Group
FD 5.58961 a
H152 5.26665 a 
summary(models$models$fit_mainEffect_aov)
##              Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo        1  0.670  0.6705   2.598 0.1171  
## PortaEnxerto  1  0.939  0.9387   3.638 0.0658 .
## Copa          2  0.133  0.0666   0.258 0.7742  
## Residuals    31  7.999  0.2580                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.44202 0.18931 28.74624 0.00000
ManejoEcologico 0.27294 0.16933 1.61194 0.11711
PortaEnxertoH152 -0.32296 0.16933 -1.90730 0.06578
CopaIAC10 -0.05720 0.20738 -0.27582 0.78452
CopaIAC5 0.09053 0.20738 0.43655 0.66546 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

BrixCorr

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "BrixCorr")
models <- fit_qf_ccsm(dframe_qf_ccsm,'BrixCorr')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 9.03333 0.15275 9.03 ± 0.15
Convencional BRF_PF+H152 9.03333 0.15275 9.03 ± 0.15
Convencional IAC10+FD 9.06667 0.15275 9.07 ± 0.15
Convencional IAC10+H152 9.00000 0.30000 9 ± 0.3
Convencional IAC5+FD 9.26667 0.05774 9.27 ± 0.06
Convencional IAC5+H152 9.30000 0.17321 9.3 ± 0.17
Ecologico BRF_PF+FD 9.26667 0.15275 9.27 ± 0.15
Ecologico BRF_PF+H152 9.33333 0.15275 9.33 ± 0.15
Ecologico IAC10+FD 9.13333 0.15275 9.13 ± 0.15
Ecologico IAC10+H152 8.86667 0.05774 8.87 ± 0.06
Ecologico IAC5+FD 9.40000 0.10000 9.4 ± 0.1
Ecologico IAC5+H152 8.96667 0.15275 8.97 ± 0.15
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Ecologico FD IAC5 9.40000 a
Ecologico H152 BRF_PF 9.33333 a
Convencional H152 IAC5 9.30000 ab
Convencional FD IAC5 9.26667 ab
Ecologico FD BRF_PF 9.26667 ab
Ecologico FD IAC10 9.13333 ab
Convencional FD IAC10 9.06667 ab
Convencional FD BRF_PF 9.03333 ab
Convencional H152 BRF_PF 9.03333 ab
Convencional H152 IAC10 9.00000 ab
Ecologico H152 IAC5 8.96667 ab
Ecologico H152 IAC10 8.86667 b 
summary(models$models$fit_full_aov)
##                          Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo                    1 0.0178 0.01778   0.711 0.40740   
## PortaEnxerto              1 0.1111 0.11111   4.444 0.04564 * 
## Copa                      2 0.2956 0.14778   5.911 0.00818 **
## Manejo:PortaEnxerto       1 0.0900 0.09000   3.600 0.06987 . 
## Manejo:Copa               2 0.2289 0.11444   4.578 0.02069 * 
## PortaEnxerto:Copa         2 0.0956 0.04778   1.911 0.16976   
## Manejo:PortaEnxerto:Copa  2 0.1067 0.05333   2.133 0.14036   
## Residuals                24 0.6000 0.02500                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 9.03333 0.09129 98.95521 0.00000
ManejoEcologico 0.23333 0.12910 1.80739 0.08325
PortaEnxertoH152 0.00000 0.12910 0.00000 1.00000
CopaIAC10 0.03333 0.12910 0.25820 0.79845
CopaIAC5 0.23333 0.12910 1.80739 0.08325
ManejoEcologico:PortaEnxertoH152 0.06667 0.18257 0.36515 0.71820
ManejoEcologico:CopaIAC10 -0.16667 0.18257 -0.91287 0.37039
ManejoEcologico:CopaIAC5 -0.10000 0.18257 -0.54772 0.58894
PortaEnxertoH152:CopaIAC10 -0.06667 0.18257 -0.36515 0.71820
PortaEnxertoH152:CopaIAC5 0.03333 0.18257 0.18257 0.85667
ManejoEcologico:PortaEnxertoH152:CopaIAC10 -0.26667 0.25820 -1.03280 0.31200
ManejoEcologico:PortaEnxertoH152:CopaIAC5 -0.53333 0.25820 -2.06559 0.04983 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico BRF_PF 9.30000 a
Convencional IAC5 9.28333 ab
Ecologico IAC5 9.18333 ab
Convencional BRF_PF 9.03333 ab
Convencional IAC10 9.03333 ab
Ecologico IAC10 9.00000 b
[[2]]
Manejo PortaEnxerto Mean Group
Ecologico FD 9.26667 a
Convencional FD 9.12222 a
Convencional H152 9.11111 a
Ecologico H152 9.05556 a
[[3]]
Copa PortaEnxerto Mean Group
IAC5 FD 9.33333 a
BRF_PF H152 9.18333 ab
BRF_PF FD 9.15000 ab
IAC5 H152 9.13333 ab
IAC10 FD 9.10000 ab
IAC10 H152 8.93333 b 
summary(models$models$fit_noInt3_aov)
##                     Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo               1 0.0178 0.01778   0.654 0.4260  
## Copa                 2 0.2956 0.14778   5.437 0.0106 *
## PortaEnxerto         1 0.1111 0.11111   4.088 0.0536 .
## Manejo:Copa          2 0.2289 0.11444   4.211 0.0261 *
## Manejo:PortaEnxerto  1 0.0900 0.09000   3.311 0.0803 .
## Copa:PortaEnxerto    2 0.0956 0.04778   1.758 0.1923  
## Residuals           26 0.7067 0.02718                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 8.96667 0.08689 103.19571 0.00000
ManejoEcologico 0.36667 0.10991 3.33612 0.00257
CopaIAC10 0.10000 0.11658 0.85782 0.39883
CopaIAC5 0.36667 0.11658 3.14533 0.00412
PortaEnxertoH152 0.13333 0.10991 1.21314 0.23598
ManejoEcologico:CopaIAC10 -0.30000 0.13461 -2.22867 0.03469
ManejoEcologico:CopaIAC5 -0.36667 0.13461 -2.72393 0.01138
ManejoEcologico:PortaEnxertoH152 -0.20000 0.10991 -1.81970 0.08034
CopaIAC10:PortaEnxertoH152 -0.20000 0.13461 -1.48578 0.14936
CopaIAC5:PortaEnxertoH152 -0.23333 0.13461 -1.73341 0.09487 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico BRF_PF 9.30000 a
Convencional IAC5 9.28333 ab
Ecologico IAC5 9.18333 ab
Convencional BRF_PF 9.03333 ab
Convencional IAC10 9.03333 ab
Ecologico IAC10 9.00000 b
[[2]]
Manejo Mean Group
Ecologico 9.16111 a
Convencional 9.11667 a
[[3]]
Copa Mean Group
IAC5 9.23333 a
BRF_PF 9.16667 ab
IAC10 9.01667 b
[[4]]
PortaEnxerto Mean Group
FD 9.19444 a
H152 9.08333 a 
summary(models$models$fit_CopPort_aov)
##                   Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo             1 0.0178 0.01778   0.503 0.4840  
## PortaEnxerto       1 0.1111 0.11111   3.142 0.0868 .
## Copa               2 0.2956 0.14778   4.179 0.0254 *
## PortaEnxerto:Copa  2 0.0956 0.04778   1.351 0.2748  
## Residuals         29 1.0256 0.03536                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 9.12778 0.08292 110.07449 0.00000
ManejoEcologico 0.04444 0.06268 0.70902 0.48397
PortaEnxertoH152 0.03333 0.10857 0.30701 0.76103
CopaIAC10 -0.05000 0.10857 -0.46052 0.64858
CopaIAC5 0.18333 0.10857 1.68858 0.10203
PortaEnxertoH152:CopaIAC10 -0.20000 0.15354 -1.30255 0.20298
PortaEnxertoH152:CopaIAC5 -0.23333 0.15354 -1.51964 0.13943 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Ecologico 9.16111 a
Convencional 9.11667 a
[[2]]
Copa Mean Group
IAC5 9.23333 a
BRF_PF 9.16667 ab
IAC10 9.01667 b
[[3]]
PortaEnxerto Mean Group
FD 9.19444 a
H152 9.08333 a 
summary(models$models$fit_mainEffect_aov)
##              Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo        1 0.0178 0.01778   0.492 0.4885  
## PortaEnxerto  1 0.1111 0.11111   3.072 0.0895 .
## Copa          2 0.2956 0.14778   4.086 0.0266 *
## Residuals    31 1.1211 0.03616                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 9.20000 0.07087 129.81073 0.00000
ManejoEcologico 0.04444 0.06339 0.70112 0.48846
PortaEnxertoH152 -0.11111 0.06339 -1.75281 0.08952
CopaIAC10 -0.15000 0.07764 -1.93207 0.06253
CopaIAC5 0.06667 0.07764 0.85870 0.39710 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

Ratio

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Ratio")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Ratio')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 1.56667 0.05774 1.57 ± 0.06
Convencional BRF_PF+H152 1.96667 0.45092 1.97 ± 0.45
Convencional IAC10+FD 1.56667 0.05774 1.57 ± 0.06
Convencional IAC10+H152 1.93333 0.05774 1.93 ± 0.06
Convencional IAC5+FD 1.66667 0.15275 1.67 ± 0.15
Convencional IAC5+H152 1.76667 0.11547 1.77 ± 0.12
Ecologico BRF_PF+FD 1.73333 0.05774 1.73 ± 0.06
Ecologico BRF_PF+H152 1.63333 0.11547 1.63 ± 0.12
Ecologico IAC10+FD 1.70000 0.26458 1.7 ± 0.26
Ecologico IAC10+H152 1.56667 0.05774 1.57 ± 0.06
Ecologico IAC5+FD 1.66667 0.05774 1.67 ± 0.06
Ecologico IAC5+H152 1.56667 0.05774 1.57 ± 0.06
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Convencional H152 BRF_PF 1.96667 a
Convencional H152 IAC10 1.93333 a
Convencional H152 IAC5 1.76667 a
Ecologico FD BRF_PF 1.73333 a
Ecologico FD IAC10 1.70000 a
Convencional FD IAC5 1.66667 a
Ecologico FD IAC5 1.66667 a
Ecologico H152 BRF_PF 1.63333 a
Convencional FD BRF_PF 1.56667 a
Convencional FD IAC10 1.56667 a
Ecologico H152 IAC10 1.56667 a
Ecologico H152 IAC5 1.56667 a 
summary(models$models$fit_full_aov)
##                          Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo                    1 0.0900  0.0900   3.115 0.09028 . 
## PortaEnxerto              1 0.0711  0.0711   2.462 0.12976   
## Copa                      2 0.0206  0.0103   0.356 0.70427   
## Manejo:PortaEnxerto       1 0.3600  0.3600  12.462 0.00171 **
## Manejo:Copa               2 0.0017  0.0008   0.029 0.97160   
## PortaEnxerto:Copa         2 0.0372  0.0186   0.644 0.53391   
## Manejo:PortaEnxerto:Copa  2 0.0450  0.0225   0.779 0.47019   
## Residuals                24 0.6933  0.0289                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.56667 0.09813 15.96511 0.00000
ManejoEcologico 0.16667 0.13878 1.20096 0.24149
PortaEnxertoH152 0.40000 0.13878 2.88231 0.00819
CopaIAC10 0.00000 0.13878 0.00000 1.00000
CopaIAC5 0.10000 0.13878 0.72058 0.47813
ManejoEcologico:PortaEnxertoH152 -0.50000 0.19626 -2.54762 0.01767
ManejoEcologico:CopaIAC10 -0.03333 0.19626 -0.16984 0.86656
ManejoEcologico:CopaIAC5 -0.16667 0.19626 -0.84921 0.40416
PortaEnxertoH152:CopaIAC10 -0.03333 0.19626 -0.16984 0.86656
PortaEnxertoH152:CopaIAC5 -0.30000 0.19626 -1.52857 0.13944
ManejoEcologico:PortaEnxertoH152:CopaIAC10 0.00000 0.27756 0.00000 1.00000
ManejoEcologico:PortaEnxertoH152:CopaIAC5 0.30000 0.27756 1.08087 0.29050 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Convencional BRF_PF 1.76667 a
Convencional IAC10 1.75000 a
Convencional IAC5 1.71667 a
Ecologico BRF_PF 1.68333 a
Ecologico IAC10 1.63333 a
Ecologico IAC5 1.61667 a
[[2]]
Manejo PortaEnxerto Mean Group
Convencional H152 1.88889 a
Ecologico FD 1.70000 ab
Convencional FD 1.60000 b
Ecologico H152 1.58889 b
[[3]]
Copa PortaEnxerto Mean Group
BRF_PF H152 1.80000 a
IAC10 H152 1.75000 a
IAC5 FD 1.66667 a
IAC5 H152 1.66667 a
BRF_PF FD 1.65000 a
IAC10 FD 1.63333 a 
summary(models$models$fit_noInt3_aov)
##                     Df Sum Sq Mean Sq F value  Pr(>F)   
## Manejo               1 0.0900  0.0900   3.169 0.08673 . 
## Copa                 2 0.0206  0.0103   0.362 0.69979   
## PortaEnxerto         1 0.0711  0.0711   2.504 0.12564   
## Manejo:Copa          2 0.0017  0.0008   0.029 0.97111   
## Manejo:PortaEnxerto  1 0.3600  0.3600  12.677 0.00145 **
## Copa:PortaEnxerto    2 0.0372  0.0186   0.655 0.52761   
## Residuals           26 0.7383  0.0284                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.59167 0.08882 17.92106 0.00000
ManejoEcologico 0.11667 0.11234 1.03848 0.30861
CopaIAC10 0.00000 0.11916 0.00000 1.00000
CopaIAC5 0.02500 0.11916 0.20980 0.83546
PortaEnxertoH152 0.35000 0.11234 3.11544 0.00444
ManejoEcologico:CopaIAC10 -0.03333 0.13759 -0.24226 0.81048
ManejoEcologico:CopaIAC5 -0.01667 0.13759 -0.12113 0.90452
ManejoEcologico:PortaEnxertoH152 -0.40000 0.11234 -3.56051 0.00145
CopaIAC10:PortaEnxertoH152 -0.03333 0.13759 -0.24226 0.81048
CopaIAC5:PortaEnxertoH152 -0.15000 0.13759 -1.09018 0.28564 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Convencional BRF_PF 1.76667 a
Convencional IAC10 1.75000 a
Convencional IAC5 1.71667 a
Ecologico BRF_PF 1.68333 a
Ecologico IAC10 1.63333 a
Ecologico IAC5 1.61667 a
[[2]]
Manejo Mean Group
Convencional 1.74444 a
Ecologico 1.64444 a
[[3]]
Copa Mean Group
BRF_PF 1.72500 a
IAC10 1.69167 a
IAC5 1.66667 a
[[4]]
PortaEnxerto Mean Group
H152 1.73889 a
FD 1.65000 a 
summary(models$models$fit_CopPort_aov)
##                   Df Sum Sq Mean Sq F value Pr(>F)
## Manejo             1 0.0900 0.09000   2.373  0.134
## PortaEnxerto       1 0.0711 0.07111   1.875  0.181
## Copa               2 0.0206 0.01028   0.271  0.765
## PortaEnxerto:Copa  2 0.0372 0.01861   0.491  0.617
## Residuals         29 1.1000 0.03793
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.70000 0.08588 19.79492 0.00000
ManejoEcologico -0.10000 0.06492 -1.54037 0.13431
PortaEnxertoH152 0.15000 0.11244 1.33400 0.19259
CopaIAC10 -0.01667 0.11244 -0.14822 0.88319
CopaIAC5 0.01667 0.11244 0.14822 0.88319
PortaEnxertoH152:CopaIAC10 -0.03333 0.15902 -0.20962 0.83543
PortaEnxertoH152:CopaIAC5 -0.15000 0.15902 -0.94328 0.35333 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Convencional 1.74444 a
Ecologico 1.64444 a
[[2]]
Copa Mean Group
BRF_PF 1.72500 a
IAC10 1.69167 a
IAC5 1.66667 a
[[3]]
PortaEnxerto Mean Group
H152 1.73889 a
FD 1.65000 a 
summary(models$models$fit_mainEffect_aov)
##              Df Sum Sq Mean Sq F value Pr(>F)
## Manejo        1 0.0900 0.09000   2.453  0.127
## PortaEnxerto  1 0.0711 0.07111   1.938  0.174
## Copa          2 0.0206 0.01028   0.280  0.758
## Residuals    31 1.1372 0.03668
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.73056 0.07138 24.24432 0.00000
ManejoEcologico -0.10000 0.06384 -1.56632 0.12743
PortaEnxertoH152 0.08889 0.06384 1.39228 0.17375
CopaIAC10 -0.03333 0.07819 -0.42630 0.67284
CopaIAC5 -0.05833 0.07819 -0.74602 0.46128 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

IT

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "IT")
models <- fit_qf_ccsm(dframe_qf_ccsm,'IT')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 1.84667 0.06807 1.85 ± 0.07
Convencional BRF_PF+H152 1.80333 0.03215 1.8 ± 0.03
Convencional IAC10+FD 1.79000 0.02646 1.79 ± 0.03
Convencional IAC10+H152 1.73333 0.04726 1.73 ± 0.05
Convencional IAC5+FD 1.84000 0.05292 1.84 ± 0.05
Convencional IAC5+H152 1.92000 0.03606 1.92 ± 0.04
Ecologico BRF_PF+FD 1.94667 0.07506 1.95 ± 0.08
Ecologico BRF_PF+H152 1.94333 0.01528 1.94 ± 0.02
Ecologico IAC10+FD 1.83667 0.03215 1.84 ± 0.03
Ecologico IAC10+H152 1.64000 0.12124 1.64 ± 0.12
Ecologico IAC5+FD 1.95333 0.04163 1.95 ± 0.04
Ecologico IAC5+H152 1.79667 0.09018 1.8 ± 0.09
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Ecologico FD IAC5 1.95333 a
Ecologico FD BRF_PF 1.94667 a
Ecologico H152 BRF_PF 1.94333 a
Convencional H152 IAC5 1.92000 a
Convencional FD BRF_PF 1.84667 ab
Convencional FD IAC5 1.84000 ab
Ecologico FD IAC10 1.83667 ab
Convencional H152 BRF_PF 1.80333 abc
Ecologico H152 IAC5 1.79667 abc
Convencional FD IAC10 1.79000 abc
Convencional H152 IAC10 1.73333 bc
Ecologico H152 IAC10 1.64000 c 
summary(models$models$fit_full_aov)
##                          Df  Sum Sq Mean Sq F value   Pr(>F)    
## Manejo                    1 0.00840 0.00840   2.283  0.14385    
## PortaEnxerto              1 0.03547 0.03547   9.637  0.00484 ** 
## Copa                      2 0.13815 0.06908  18.768 1.24e-05 ***
## Manejo:PortaEnxerto       1 0.02834 0.02834   7.699  0.01053 *  
## Manejo:Copa               2 0.03651 0.01825   4.959  0.01575 *  
## PortaEnxerto:Copa         2 0.01871 0.00935   2.541  0.09977 .  
## Manejo:PortaEnxerto:Copa  2 0.02957 0.01479   4.017  0.03127 *  
## Residuals                24 0.08833 0.00368                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.84667 0.03503 52.72207 0.00000
ManejoEcologico 0.10000 0.04953 2.01878 0.05482
PortaEnxertoH152 -0.04333 0.04953 -0.87480 0.39035
CopaIAC10 -0.05667 0.04953 -1.14398 0.26392
CopaIAC5 -0.00667 0.04953 -0.13459 0.89406
ManejoEcologico:PortaEnxertoH152 0.04000 0.07005 0.57100 0.57331
ManejoEcologico:CopaIAC10 -0.05333 0.07005 -0.76133 0.45388
ManejoEcologico:CopaIAC5 0.01333 0.07005 0.19033 0.85065
PortaEnxertoH152:CopaIAC10 -0.01333 0.07005 -0.19033 0.85065
PortaEnxertoH152:CopaIAC5 0.12333 0.07005 1.76057 0.09105
ManejoEcologico:PortaEnxertoH152:CopaIAC10 -0.18000 0.09907 -1.81690 0.08174
ManejoEcologico:PortaEnxertoH152:CopaIAC5 -0.27667 0.09907 -2.79265 0.01010 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico BRF_PF 1.94500 a
Convencional IAC5 1.88000 ab
Ecologico IAC5 1.87500 ab
Convencional BRF_PF 1.82500 bc
Convencional IAC10 1.76167 bc
Ecologico IAC10 1.73833 c
[[2]]
Manejo PortaEnxerto Mean Group
Ecologico FD 1.91222 a
Convencional FD 1.82556 ab
Convencional H152 1.81889 b
Ecologico H152 1.79333 b
[[3]]
Copa PortaEnxerto Mean Group
BRF_PF FD 1.89667 a
IAC5 FD 1.89667 a
BRF_PF H152 1.87333 a
IAC5 H152 1.85833 a
IAC10 FD 1.81333 a
IAC10 H152 1.68667 b 
summary(models$models$fit_noInt3_aov)
##                     Df  Sum Sq Mean Sq F value   Pr(>F)    
## Manejo               1 0.00840 0.00840   1.853  0.18512    
## Copa                 2 0.13815 0.06908  15.232 4.18e-05 ***
## PortaEnxerto         1 0.03547 0.03547   7.822  0.00958 ** 
## Manejo:Copa          2 0.03651 0.01825   4.025  0.03000 *  
## Manejo:PortaEnxerto  1 0.02834 0.02834   6.249  0.01907 *  
## Copa:PortaEnxerto    2 0.01871 0.00935   2.062  0.14744    
## Residuals           26 0.11791 0.00453                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.80861 0.03549 50.95842 0.00000
ManejoEcologico 0.17611 0.04489 3.92281 0.00057
CopaIAC10 -0.01167 0.04762 -0.24501 0.80837
CopaIAC5 0.06250 0.04762 1.31255 0.20080
PortaEnxertoH152 0.03278 0.04489 0.73011 0.47185
ManejoEcologico:CopaIAC10 -0.14333 0.05498 -2.60683 0.01493
ManejoEcologico:CopaIAC5 -0.12500 0.05498 -2.27340 0.03150
ManejoEcologico:PortaEnxertoH152 -0.11222 0.04489 -2.49971 0.01907
CopaIAC10:PortaEnxertoH152 -0.10333 0.05498 -1.87934 0.07145
CopaIAC5:PortaEnxertoH152 -0.01500 0.05498 -0.27281 0.78716 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Ecologico BRF_PF 1.94500 a
Convencional IAC5 1.88000 a
Ecologico IAC5 1.87500 a
Convencional BRF_PF 1.82500 a
Convencional IAC10 1.76167 ab
Ecologico IAC10 1.73833 b
[[2]]
Manejo Mean Group
Ecologico 1.85278 a
Convencional 1.82222 a
[[3]]
Copa Mean Group
BRF_PF 1.8850 a
IAC5 1.8775 a
IAC10 1.7500 b
[[4]]
PortaEnxerto Mean Group
FD 1.86889 a
H152 1.80611 b 
summary(models$models$fit_CopPort_aov)
##                   Df  Sum Sq Mean Sq F value   Pr(>F)    
## Manejo             1 0.00840 0.00840   1.333 0.257623    
## PortaEnxerto       1 0.03547 0.03547   5.629 0.024519 *  
## Copa               2 0.13815 0.06908  10.961 0.000285 ***
## PortaEnxerto:Copa  2 0.01871 0.00935   1.484 0.243402    
## Residuals         29 0.18275 0.00630                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.88139 0.03500 53.74702 0.00000
ManejoEcologico 0.03056 0.02646 1.15474 0.25762
PortaEnxertoH152 -0.02333 0.04583 -0.50911 0.61453
CopaIAC10 -0.08333 0.04583 -1.81825 0.07937
CopaIAC5 0.00000 0.04583 0.00000 1.00000
PortaEnxertoH152:CopaIAC10 -0.10333 0.06482 -1.59426 0.12172
PortaEnxertoH152:CopaIAC5 -0.01500 0.06482 -0.23143 0.81861 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Ecologico 1.85278 a
Convencional 1.82222 a
[[2]]
Copa Mean Group
BRF_PF 1.8850 a
IAC5 1.8775 a
IAC10 1.7500 b
[[3]]
PortaEnxerto Mean Group
FD 1.86889 a
H152 1.80611 b 
summary(models$models$fit_mainEffect_aov)
##              Df  Sum Sq Mean Sq F value   Pr(>F)    
## Manejo        1 0.00840 0.00840   1.293 0.264202    
## PortaEnxerto  1 0.03547 0.03547   5.458 0.026119 *  
## Copa          2 0.13815 0.06908  10.629 0.000305 ***
## Residuals    31 0.20145 0.00650                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.90111 0.03004 63.28022 0.00000
ManejoEcologico 0.03056 0.02687 1.13712 0.26420
PortaEnxertoH152 -0.06278 0.02687 -2.33626 0.02612
CopaIAC10 -0.13500 0.03291 -4.10208 0.00027
CopaIAC5 -0.00750 0.03291 -0.22789 0.82123 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))

VitC

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "VitC")
models <- fit_qf_ccsm(dframe_qf_ccsm,'VitC')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 34.15667 2.19805 34.16 ± 2.2
Convencional BRF_PF+H152 37.52333 1.58418 37.52 ± 1.58
Convencional IAC10+FD 34.25333 1.33121 34.25 ± 1.33
Convencional IAC10+H152 33.71000 2.57874 33.71 ± 2.58
Convencional IAC5+FD 32.83667 2.87114 32.84 ± 2.87
Convencional IAC5+H152 33.98333 1.23678 33.98 ± 1.24
Ecologico BRF_PF+FD 31.96000 1.10585 31.96 ± 1.11
Ecologico BRF_PF+H152 33.32333 3.61777 33.32 ± 3.62
Ecologico IAC10+FD 34.23000 2.17690 34.23 ± 2.18
Ecologico IAC10+H152 36.66667 0.76846 36.67 ± 0.77
Ecologico IAC5+FD 31.53667 1.24340 31.54 ± 1.24
Ecologico IAC5+H152 33.23333 1.44901 33.23 ± 1.45
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
Manejo PortaEnxerto Copa Mean Group
Convencional H152 BRF_PF 37.52333 a
Ecologico H152 IAC10 36.66667 ab
Convencional FD IAC10 34.25333 ab
Ecologico FD IAC10 34.23000 ab
Convencional FD BRF_PF 34.15667 ab
Convencional H152 IAC5 33.98333 ab
Convencional H152 IAC10 33.71000 ab
Ecologico H152 BRF_PF 33.32333 ab
Ecologico H152 IAC5 33.23333 ab
Convencional FD IAC5 32.83667 ab
Ecologico FD BRF_PF 31.96000 ab
Ecologico FD IAC5 31.53667 b 
summary(models$models$fit_full_aov)
##                          Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo                    1   7.60   7.599   1.868 0.1844  
## PortaEnxerto              1  22.40  22.404   5.507 0.0275 *
## Copa                      2  21.33  10.665   2.621 0.0934 .
## Manejo:PortaEnxerto       1   0.58   0.583   0.143 0.7084  
## Manejo:Copa               2  32.69  16.347   4.018 0.0313 *
## PortaEnxerto:Copa         2   3.13   1.564   0.384 0.6850  
## Manejo:PortaEnxerto:Copa  2   9.31   4.657   1.145 0.3351  
## Residuals                24  97.64   4.069                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 34.15667 1.16455 29.33034 0.00000
ManejoEcologico -2.19667 1.64692 -1.33380 0.19479
PortaEnxertoH152 3.36667 1.64692 2.04422 0.05205
CopaIAC10 0.09667 1.64692 0.05870 0.95368
CopaIAC5 -1.32000 1.64692 -0.80149 0.43071
ManejoEcologico:PortaEnxertoH152 -2.00333 2.32910 -0.86013 0.39823
ManejoEcologico:CopaIAC10 2.17333 2.32910 0.93312 0.36005
ManejoEcologico:CopaIAC5 0.89667 2.32910 0.38498 0.70364
PortaEnxertoH152:CopaIAC10 -3.91000 2.32910 -1.67876 0.10617
PortaEnxertoH152:CopaIAC5 -2.22000 2.32910 -0.95316 0.35001
ManejoEcologico:PortaEnxertoH152:CopaIAC10 4.98333 3.29385 1.51292 0.14336
ManejoEcologico:PortaEnxertoH152:CopaIAC5 2.55333 3.29385 0.77518 0.44580 
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Convencional BRF_PF 35.84000 a
Ecologico IAC10 35.44833 a
Convencional IAC10 33.98167 a
Convencional IAC5 33.41000 a
Ecologico BRF_PF 32.64167 a
Ecologico IAC5 32.38500 a
[[2]]
Manejo PortaEnxerto Mean Group
Convencional H152 35.07222 a
Ecologico H152 34.40778 a
Convencional FD 33.74889 a
Ecologico FD 32.57556 a
[[3]]
Copa PortaEnxerto Mean Group
BRF_PF H152 35.42333 a
IAC10 H152 35.18833 a
IAC10 FD 34.24167 a
IAC5 H152 33.60833 a
BRF_PF FD 33.05833 a
IAC5 FD 32.18667 a 
summary(models$models$fit_noInt3_aov)
##                     Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo               1   7.60   7.599   1.847 0.1858  
## Copa                 2  21.33  10.665   2.593 0.0940 .
## PortaEnxerto         1  22.40  22.404   5.446 0.0276 *
## Manejo:Copa          2  32.69  16.347   3.974 0.0312 *
## Manejo:PortaEnxerto  1   0.58   0.583   0.142 0.7097  
## Copa:PortaEnxerto    2   3.13   1.564   0.380 0.6875  
## Residuals           26 106.96   4.114                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 34.78472 1.06898 32.53997 0.00000
ManejoEcologico -3.45278 1.35217 -2.55351 0.01687
CopaIAC10 -1.14917 1.43419 -0.80126 0.43024
CopaIAC5 -1.95833 1.43419 -1.36546 0.18381
PortaEnxertoH152 2.11056 1.35217 1.56087 0.13065
ManejoEcologico:CopaIAC10 4.66500 1.65606 2.81692 0.00914
ManejoEcologico:CopaIAC5 2.17333 1.65606 1.31235 0.20087
ManejoEcologico:PortaEnxertoH152 0.50889 1.35217 0.37635 0.70971
CopaIAC10:PortaEnxertoH152 -1.41833 1.65606 -0.85645 0.39957
CopaIAC5:PortaEnxertoH152 -0.94333 1.65606 -0.56962 0.57382 
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
Manejo Copa Mean Group
Convencional BRF_PF 35.84000 a
Ecologico IAC10 35.44833 a
Convencional IAC10 33.98167 a
Convencional IAC5 33.41000 a
Ecologico BRF_PF 32.64167 a
Ecologico IAC5 32.38500 a
[[2]]
Manejo Mean Group
Convencional 34.41056 a
Ecologico 33.49167 a
[[3]]
Copa Mean Group
IAC10 34.71500 a
BRF_PF 34.24083 a
IAC5 32.89750 a
[[4]]
PortaEnxerto Mean Group
H152 34.74000 a
FD 33.16222 b 
summary(models$models$fit_CopPort_aov)
##                   Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo             1   7.60   7.599   1.571 0.2200  
## PortaEnxerto       1  22.40  22.404   4.633 0.0398 *
## Copa               2  21.33  10.665   2.206 0.1283  
## PortaEnxerto:Copa  2   3.13   1.564   0.323 0.7263  
## Residuals         29 140.24   4.836                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 33.51778 0.96968 34.56581 0.00000
ManejoEcologico -0.91889 0.73301 -1.25358 0.22001
PortaEnxertoH152 2.36500 1.26961 1.86278 0.07265
CopaIAC10 1.18333 1.26961 0.93205 0.35901
CopaIAC5 -0.87167 1.26961 -0.68656 0.49781
PortaEnxertoH152:CopaIAC10 -1.41833 1.79550 -0.78994 0.43598
PortaEnxertoH152:CopaIAC5 -0.94333 1.79550 -0.52539 0.60331 
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
Manejo Mean Group
Convencional 34.41056 a
Ecologico 33.49167 a
[[2]]
Copa Mean Group
IAC10 34.71500 a
BRF_PF 34.24083 a
IAC5 32.89750 a
[[3]]
PortaEnxerto Mean Group
H152 34.74000 a
FD 33.16222 b 
summary(models$models$fit_mainEffect_aov)
##              Df Sum Sq Mean Sq F value Pr(>F)  
## Manejo        1   7.60   7.599   1.643 0.2094  
## PortaEnxerto  1  22.40  22.404   4.845 0.0353 *
## Copa          2  21.33  10.665   2.306 0.1165  
## Residuals    31 143.36   4.625                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 33.91139 0.80144 42.31301 0.00000
ManejoEcologico -0.91889 0.71683 -1.28188 0.20939
PortaEnxertoH152 1.57778 0.71683 2.20105 0.03531
CopaIAC10 0.47417 0.87793 0.54009 0.59299
CopaIAC5 -1.34333 0.87793 -1.53011 0.13613 
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))